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ABSTRACT 

Using Spitzer IRAC and MIPS observations of the Large Magellanic Cloud, we have identified 13 
objects that have extremely red mid-IR colors. Follow-up Spitzer IRS observations of seven of these 
sources reveal varying amounts of SiC and C2H2 absorption as well as the presence of a broad MgS 
feature in at least two cases, indicating that these are extreme carbon stars. Preliminary estimates 
find these objects have luminosities of 4-11x10'^ L© and preliminary model fitting gives mass-loss 
rates between 4xl0~^ and 2x10""* Mq yr~^, higher than any known carbon-rich AGB star in the 
LMC. These spectral and physical properties require careful reconsideration of dust condensation and 
mass-loss processes for carbon stars in low metallicity environments. 

Subject headings: Magellanic Clouds — stars: AGB and post- AGB — stars: carbon — infrared: stars 



1. INTRODUCTION 

We have undertaken a study of star formation in th e 
Large Magellanic Cloud fLMC: ICniendl fc Chul[200l . 
using archival Spitzer Space Telescope observations, such 
as those of Su rveying the Agents of a Galaxy's Evo- 
lution fSAGE: [Meixner et al.ll2006D . In the process of 
identifying young stellar objects (YSOs) we noticed 13 
bright mid-IR sources (Table 1) that all had similar pho- 
tometric properties with spectral energy distributions 
(SEDs) markedly different from those of evolved stars 
or YSOs: (1) they have extremely red mid-IR colors, 
[4.5]-[8.0]>4.0; (2) their SEDs, peaking between 8 and 
24 /im, can be moderately well fit by blackbodies with 
effective temperatures of ~230-320 K; (3) they all fall 
in a narrow range of brightness, 7.0>[8.0]>8.5; and (4) 
none of the sources have counterparts in the Two Mi- 
cron All Sky Survey P oint Source Catalog (2MASS PSC; 
iSkrutskie etalMm or in the Digitized Sky Survey. We 
dubbed these sources Extremely Red Objects (EROs). 

The EROs are not likely asteroids because of their 
high brightnesses and lack of proper motions in 
SAGE observations from two epochs separated by '^3 
months. They cannot be background galaxies or Galac- 
tic sources, as t he Spitzer Wide-area IR Extragalactic 
burvey (SWIRE; iLonsdale et"alll2003f ) and the Galactic 
Legacy IR Mid-Plane Survey Extraordinaire (GLIMPSE; 
iBeni amin et al]l2003f ) do not have counterparts with sim- 
ilar mid-IR colors and brightnesses. These EROs are 
most likely associated with the LMC, and their high lu- 
minosities imply that their Galactic counterparts would 
have saturated in the SWIRE and GLIMPSE Surveys. 

Queries of the SIMBAD database found that some 
EROs had been detected previously in either MSX or 
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IRAS observations (e.g., iSchwerin"^ 119890 . Five have 
been suggested as possible obscured AGB stars by 
iLoup et al] (|1997D : however, many of our EROs have also 
been suggested to be YSO candidates (jWhitnev et al.l 
120081 ). As these EROs have never been confirmed or 
rejected as YSOs spectroscopically, we included seven 
among our follow-up observations of massive YSOs in 
the LMC using the Spitzer InfraRed Spectrograph (IRS; 
iHouck et al.l |2004 ). These IRS observations show un- 
ambiguously that the EROs are carbon stars; further- 
more, the spectra reveal silicon carbide (SiC) absorption 
features. While SiC emission features are fairly com- 
mon in both Galactic and LMC carbon stars, this is the 
first clear detection of SiC absorption for LMC carbon 
stars. Preliminary analysis of the IRS spectra suggests 
that these are extraordinary carbon stars with very high 
mass-loss rates. This paper reports their discovery. In §2 
we describe their basic photometric properties and their 
possible optical and near-IR counterparts. In §3 we in- 
troduce our IRS observations and in §4 we discuss the 
results. 

2. PHOTOMETRIC OBSERVATIONS 

2.1. Mid- Infrared Observations 

Spitzer observations of the EROs, from SAGE 
(Meixner et al. 2006) . were made with the InfraRed Ar- 
ray Camera (IRAC; iPazio et al] [200 41 and t he Multi- 
band Imaging Photometer for Spitzer (MIPS: ,Rieke et al.l 
|2004[ ). Aperture photometry of the EROs was obtained 
for the IRAC 3.6, 4.5, 5.8, and 8.0 /im bands and 
the MIPS 24 and 70 /im bands using the task phot 
in IRAF. A detailed description of the data reduction, 
p hotometric extraction, a nd uncertainties can be found 
in iGruendl fc Chul (|2008f) . When multiple observations 
were available, the photometric measurements were aver- 
aged. Table [1] presents the average photometric results: 
column 1 lists source names that provide epoch 2000 co- 
ordinates, columns 2-7 give flux densities and uncertain- 
ties at 3.6, 4.5, 5.8, 8.0, 24, and 70 /xm, and column 
8 gives cross-identifications to previously known sources. 
Figure [1] plots the locations of the EROs in a mid-IR [8.0] 
vs. [4.5] — [8.0] color-magnitude diagram along with other 
red so urces classified as evolved stars bv IGruendl fc Chul 
(|2008l ). The EROs are at the extreme red end of the tail 
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TABLE 1 

Mid-Infrared Photometric Properties and Derived Parameters 



Cross 



Mass- Loss 



Source ID 




m4.5 




ms.O 


m24.0 


m70.o 


Identification 


Lbol 


Rate 








[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 






[Lq] 


[Mq yr" 


1] 


050231.49 


-680535.8 


15.92±.08 


12.75±.05 


9.91±.05 


7.26±.05 


3.06±.ll 


2.09±.23 






7800 


1.1x10- 


4 


050343.02 


-664456.7 


14.92±.06 


12.93±.06 


10.68±.06 


8.25±.06 


2.96±.ll 


1.71±.24 


IKAb 05036- 


'6649 


7000 


1.5x10" 


'4 


050405.60 


-682340.3 


16.47±.10 


13.12±.06 


10.27±.05 


7.65±.05 


3.63±.ll 


>3.32^ 


IRAS 05042- 


'6827 


4950 


6.3x10" 


-5 


051301.75 


-693351.0 


15.87±.08 


13.70±.06 


10.69±.06 


7.78±.06 


3.35±.ll 


>1.59=' 


IRAS 05133- 


'6937 


5800 


8.9x10" 


'5 


051811.70 


-703027.0 


14.97±.07 


12.26±.06 


9.89±.05 


7.61±.05 


2.82±.ll 


1.59±.22 


IRAS 05187- 


'7033 


8850 


1.3x10" 


'4 


051848.36 


-693334.7 


14.41±.07 


12.34±.06 


9.72±.05 


7.11±.05 


3.32±.ll 


1.24±.24 


IRAS 05191- 


'6936 


7750 


8.3x10" 


'5 


052540.63 


-700827.2 


16.02±.13 


13.20±.06 


10.45±.06 


7.80±.05 


3.82±.ll 


2.58±.26 


IRAS 05260- 


'7010 


4300 


5.9x10" 


'5 


052937.89 


-724952.9 


13.61±.06 


11.14±.05 


9.02±.05 


7.09±.05 


3.81±.ll 


2.92±.26 


IRAS 05305- 


'7251 


5550 


4.2x10" 


'5 


053044.10 


-714300.5 


13.47±.06 


12.13±.05 


10.21±.05 


7.76±.06 


2.75±.ll 


1.13±.22 


IRAS 05315- 


'7145 


9500 


1.7x10" 


'4 


054134.73 


-694209.3 


14.84±.07 


11.91±.05 


9.40±.05 


7.11±.05 


3.50±.ll 


>2.35^ 


IRAS 05420- 


'6943 


6400 


6.6x10" 


'5 


054859.98 


-703322.5 


15.33±.08 


13.70±.06 


11.32±.06 


8.41±.06 


2.41±.ll 


0.72±.22 


IRAS 05495- 


'7034 


11100 


2.3x10" 


'4 


055026.08 


-695603.1 


14.26±.05 


11.81±.05 


9.25±.05 


6.94±.05 


3.10±.ll 


2.36±.23 


IRAS 05509- 


'6956 


8250 


9.0x10" 


'5 


055133.60 


-711933.9 


14.06±.06 


11.28±.05 


9.03±.05 


7.01±.05 


3.87±.ll 


3.38±.26 


IRAS 05522- 


'7120 


5200 


4.4x10" 


'5 


^ For non- 


detections a 


3-(T upper-limit is given. 
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Fig. 1. — Color-magnitude diagram showing IRAC [8.0] vs. 
[4.5] — [8.0]. The gray scale Hess diagram shows all data from the 
LMC. The location of the EROs are plotted with filled circles while 
other sources classified as evolved stars by . Gruendl fc Chu (2008 ) 
have been plotted with open circles. 

formed by the evolved stars. 

2.2. Complementary Near-Infrared and Optical 
Photometry 

None of the EROs have a counterpart in the 2JVIASS 
PSC. Thus, we obtained deeper ncar-IR observations 
at J and iCc-bands for six of the EROs usine; the IR 



Side Port Imager (ISPI: lvan der BUek et all200l) on the 
Blanco 4m telescope at the Cerro Tololo Inter- American 
Observatory in 2007 February. The resulting images have 
a typical effective exposure time of ~300 s and 600 s in 
the J and Xj-bands, respectively, and are flux-calibrated 
using stars in the 2IVIASS PSC. The angular resolution of 
the images is typically < I'.'O and the astrometric accu- 
racy is better than 0'.'2. For a m ore detailed descr i ption of 
the near-IR data reduction, see iGruendl fc CkH (l2008h . 

Table [5] presents the near-IR photometry and 3-ct 
upper limits for the six EROs observed. Only two 
EROs, 051301.75-693351.0 and 053044.10-714300.5, 
have near-IR counterparts projected within 0'.'2 as indi- 



TABLE 2 

Possible Near-Infrared Counterparts 



Source ID 




mK 


Offset 






[mag] 


[mag] 


["] 


ISPI Observations 


050231.49- 


680535.8 


>19.67 


>18.63 




050405.60- 


682340.3 


>19.84 


>18.60 




051301.75- 


693351.0 


16.91±0.05 


16.65±0.06 


<0.2 


052937.89- 


724952.9 


>19.52 


>18.74 




053044.10- 


714300.5 


17.79±0.07 


16.77±0.07 


<0.2 


055026.08- 


695603.1 


>19.82 


>18.74 






Possible Counterparts in 


IRSF Catalog 




051301.75- 


693351.0 


16.98±0.02 


16.27±0.06 


0.1 


051848.36- 


693334.7 


17.27±0.03=' 


16.66±0.12=' 


1.0 


053044.10- 


714300.5 


18.54±0.08 


17.23±0.22 


0.2 


052540.63- 


700827.2 


17.77±0.04=' 


17.08±0.17'' 


0.7 


054134.73- 


694209.3 


19.61±0.21=' 


16.70±0.09'' 


0.6 



^ Doubtful counterpart 

cated in the last column of the table. To supplement our 
ISPI observations, we have used the JVIagellanic Clouds 
Point Source Catalog from the InfraRed Survey Facil- 
ity fIRSF: IKato et al.li2007t) . This survey has an astro- 
metric accuracy of ~0'.'l and 90% completeness limits of 
18.5 and 17.4 for mj and uik, respectively. IRSF sources 
with spatial coincidence of 1" or better are also included 
in Table [21 Only two IRSF sources are within 0'.'2 or 
less from their respective EROs and are the same coun- 
terparts identified by our ISPI observations. The other 
three IRSF sources are included for completeness but are 
unlikely to be near-IR counterparts to the EROs. 

To search for optical counterparts, we used the opti- 
cal photometry f rom the Magellanic C louds Photometric 
burvey (JVICPS; IZaritskv et all [200l which has an as- 
trometric accuracy of ~0'.'5 and is generally complete to 
mv =20. We find that three EROs have faint (my >19) 
optical sources within 1", but are likely unrelated as none 
have a spatial coincidence better than 0'.'7. In Figure [2] 
we combine all photometric data for each source to get 
a broad SED extending from the optical to 70 /im. 

3. IRS OBSERVATIONS 

Seven EROs were included in our Spitzer IRS survey 
of YSOs (PID=40650), which used the Short-Low (SL) 
modules to obtain spectra over '^5.2-14.5 /im with spec- 
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Fig. 2. — The SEDs of the EROs using optical photometry 
from the MCPS, near-IR photometry from our ISPI observations 
or the IRSF survey, and mid-IR photometry from IRAC and MIPS 
measurements. Measurements with good astrometric coincidence 
are plotted with filled symbols, while those with poorer coincidence 
are indicated by open symbols. Upper limits for non-detections are 
shown as arrows with a crossbar at the 3-cr significance limit. 

tral resolution R (= 6X/X) of 64-128, and the Short- 
High (SH) and Long-High (LH) modules to obtain spec- 
tra over ~9.9-37.2 fim with spectral resolution of ~600. 
Background observations with the same integration times 
were made in nearby regions with low surface bright- 
nesses in the IRAS 12 and 25 fim maps. 

To extract spectra, we used the basic calibrated 
data (BCD) from the Spitzer Science Center's pipeline. 
Rogue pixels and flagged data were cleaned using the 
IRSCLEAN package (ver 1.9). Multiple exposures for 
each slit position were median averaged. The background 
was then removed in each module by subtracting the 
cleaned BCD image data at the background position. 
For the SL modules, an additional local background was 
subtracted by differencing the two nod positio ns. Spectra 
were extracted using the SMART software (Hi gdon et al.l 
l2004f ) with a full aperture extraction for the SH and LH 
modules and the tapered column point source method 
for the SL modules. The SH and LH spectra were auto- 
matically defringed using the IRSFRINGE package (ver 
1.1). To combine all spectra, we used the SL extracted 
spectra to set the flux level (as their background sub- 
traction was more robust) and then applied a multiplica- 
tive scale factor for the SH and LH spectra based on the 
continuum where the spectra overlapped. The resulting 
spectra span the wavelength range from 5.2 to 37.2 jum. 
A mo re detailed description is presented in ' Scale et al.l 
(|2009[ ) which reports the results of our IRS survey of 
LMC YSO candidates. 

4. RESULTS 

Figure [3] presents the IRS spectra of seven EROs along 
with the IRAC and MIPS mid-IR photometry. Clearly, 
these spectra are dominated by dust continuum emis- 
sion. Interestingly, the 11.3 /im SiC feature is detected 
in absorption in nearly all spectra; in Figure [3] the spec- 
tra have been ordered from top to bottom based on the 
strength of the SiC absorption feature. In addition, the 
C2H2 (acetylene) 13.7 fiia absorption feature appears in 



T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 
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10 20 30 
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Fig. 3. — IRS spectra for seven EROs. The spectra are roughly 
ordered from top to bottom based on the increasing strength of 
the SiC feature at 11.3 fj,m. Each spectrum has been offset by a 
constant value in Jy which is given after the source name in the 
figure. Our mid-IR photometry for each source are plotted with 
each spectrum for comparison using alternating symbols. 

every spectrum, and the C2H2 7.5 /im absorption fea- 
ture appears in some spectra. In contrast, the MgS 
emission feature at ^26 fim is clearly detected in the 
spectrum of 055026.08—695603.1 and possibly in that 
of 052937.89-724952.9. Finally, we note that the only 
source without SiC absorption, 051811.70-703027.0, has 
a nearly featureless spectrum similar to AFGL618 (the 
Westbrook Nebula) a carbon-rich prot oplanctary nebula 
(jKraemer et al.ll2002t ICernicharo et al..,2001a .b). 

4.1. Mass-Loss Rates 

The extreme red colors of these EROs, [8.0] -[24]= 3- 
5, suggest that they have high opacities and hence high 
mass-loss rates, M. Simple fits to the mid-IR SEDs of 
the EROs give bolometric luminosit ies in the ra nge of 4- 
llxlO^Lg (for a distance of 50 kpc;'Feast"1999). Prclim- 
inary analysis using the radiative transfer model DUSTY 
(jlvezic fc Elitzud Il995l ) shows that the highest optical 
depths (r) are found for 054859.98-703322.5, with r's 
of 270, 8.1, and 10.0 at 1, 10 and 11.3 /im, respectively. 
Assuming a gas-to-dust ratio of 200, the estimated M are 
in the range of 0.4-2.3 x 10"'' M© yr'^ (see Tabled]). We 
have adjusted values of M from other works to reflect a 
gas-to-dust ratio of 200. 

Our derived values for Al are higher than previously 
known for carbon-rich AGB stars in the LMC. Fur- 
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thermore, these 13 sources more than double the pre- 
viously known 8 carbon-rich AGB stars in the LMC 
with M h igher t han 10~ ^ A/^ yr"-'^ ()van Loon et al.l 
119991: iGroe ncwege n et al.l 12007.) . More intriguingly, 
these sources typically exceed the maximum M ex- 
pected for both oxygen-rich and carbon-rich AGB 
stars (jvan Loon et al.l Il999). For example, the M of 
054859.98-703322.5, 2.3x10"'' Mq yr"!, is much higher 
than the maximum expected for its bolometric luminos- 
ity, 4 X 10~5 Mq yr-i. 

Comparisons of carbon stars in the Magellanic Clouds 
and t he Galaxy show that M does not depend on metal- 
hcity (jMatsuura et al.ll2005t ISloan et al.ll2008f l. because 
M depends on the amount of carbon available in the at- 
mosphere, which is ultimately synthesized within these 
stars. These new extreme carbon stars show that the 
highest M found in the LMC ar e comparable to those of 
their Galactic counterparts fe.g.. lSpeck et al.1 ^008'). The 
mean bolometric luminosity of the EROs, ~7,100 Lq, 
corresponds to main-sequence masses o f 1.5-2.5 Mq (us- 
ing Fig. 20 in lVassihadis fc Woodlll993h . where the lower 
limit is set by the mass range to become carbon stars 
(I.5-5.OM0). It is puzzling that such low-mass stars 
could reach the highest M. This might be caused by 
some specific combinations of abundance and pressure- 
temperature profile in the stellar at mosphere, resulting 
in an unusual condensation sequence (jLodders fc FeglevI 
Il995| ). It is also possible that a carbon-rich dust disk is 
responsible for the SiC absorption then the M based on 
a spherical symmetry is not valid. These possibilities will 
be considered in greater detail in future work. 

4.2. Spectral Features 

The occurrence of SiC absorption features is sur- 
prising. These features are rare in Galactic extreme 
carbon stars, with only two unambiguous detections 



(jSpeck et al.l [19971) and another 8 more tentative de- 
tectio ns (jPitman. Speck fc Hofmeisteill200dlSpeck et al.l 
[2OO8I) . This raises the question of why the SiC absorp- 
tion feature is seemingly common in the LMC extreme 
carbon stars. 

Recent Spitzer observations show that planetary neb- 
ulae (PNe) in the LMC often show strong SiC emis - 
sion (jStanghellini et al.ll2007l : iBernard-Salas et al.ll2008D. 



while in Galactic PN e, it is generally weak ( Barlowlll983l : 
iRinehart et"ani2002f) . The presence of SiC features in 
EROs and PNe in the LMC suggests that (a) SiC grains 
are formed more ab undantly in the LM C than previ- 
ous recognized (e.g.. IZijlstra et ani2006l ): and (b) that 
the LMC metallicity allows SiC to form and survive in 
the late AGB phase and into the PN phase more effi- 
cie ntly than in the Galax y. It has been hypothesized 
bv iLeisenring et aTl ()2008f ) that MgS tends to form as a 
mantle on pre-existing SiC grains and gives rise to the 
emission feature at 30 /im. This hypothesis is supported 
by a recent study of the correlation between mass-loss 
rate and the amount of SiC r equired to model G alactic 
extreme carbon star spectra (jSpeck et al.ll2008l ). Only 
two of the LMC EROs show evidence for a MgS fea- 
ture at 30 /im indicating that the SiC grains may not 
be coated for some of these objects. Current models for 
the effect of metallicity on condensation s equences in the 
Galaxy and the Magellanic Clouds (e.g., Lagadec et al.l 
120071: (Leisenring et all 120081 : ISpeck et al.ll2006D needto 
be re-analyzed in light of this discovery. 
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